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Introduction: 

A  long-standing  clinical  challenge  in  the  management  of  prostate  cancer  is  its  heterogeneous  response  to 
androgen  deprivation  therapy  (ADT),  a  standard  treatment  to  disrupt  the  androgen  receptor  (AR)  signaling 
pathway,  since  AR  has  a  profound  effect  on  prostate  cancer  development  through  the  regulation  of  not  only 
transcriptional  networks  but  also  genomic  stability,  and  the  development  of  gene  fusions[l].  While  ADT  is 
effective  in  most  patients  with  prostate  cancer,  prostate  cancers  inevitably  become  resistant  to  ADT  and 
castration-resistant  prostate  cancer  (CRPC)  emerges [2],  Since  multiple  molecular  mechanisms  modulating 
resistance  to  ADT  have  been  proposed  [3],  the  recent  surge  in  prostate  cancer  genomic  information  may  permit 
molecular  classification  of  prostate  cancer  [4].  As  shown  in  Fig.  la,  whole  exome  sequencing  of  metastatic 
castration  resistant  prostate  tumors  (mCRPC)  from  the  Michigan  cohort  [5]  and  the  East  Coast  Stand  Up  to 
Cancer/Prostate  Cancer  Foundation  (SU2C/PCF)  study  recently  revealed  that  5-10%  of  tumors  harbor 
promyelocytic  leukemia  zinc  finger  (PFZF)  focal  homozygous  deletions.  PLZF  also  known  as  the  zinc  finger 
and  BTB-containing  protein  16  (ZBTB16)  was  originally  identified  as  a  gene  fused  to  RARa  in  acute 
promyelocytic  leukemia  (APF)  patients  [6]  and  was  reported  as  an  androgen-responsive  tumor  suppressor 
gene[7].  However,  the  role  of  PFZF  in  prostate  cancer  progression  is  unclear. 

Keyword: 
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Accomplishment : 

I  interrogated  the  biological  consequences  of  the  loss  of  PLZF  expression  and  performed  a  bioinformatic 
analysis  of  PLZF  ChIP  seq  and  PLZF-regulated  gene  expression  profiling  to  define  the  PLZF  transcriptional 
program  as  shown  in  the  Fig.  4  below.  In  addition,  I  examined  the  biological  function  of  PLZF  in  CRPC 
development  in  vitro  and  in  vivo  as  shown  in  the  Fig.  1,  2  and  3.  Collectively,  our  preliminary  data  strongly 
support  that  PLZF  functions  as  a  tumor  suppressor  in  prostate  cancer  and  demonstrate  that  PLZF  inactivation  as 
an  important  molecular  event  for  CRPC  development. 

What  were  the  major  goals  of  the  project? 

Task  1.  Identify  PLZF  transcription  program  mediating  castration  resistance  in  prostate  cancer  cells. 
(Months  1-12) 

a.  Determine  the  PLZF-regulated  genes  in  LNCaP  cells  by  using  the  microarray  gene  profiling 
technique  (Months  1-4)  and  compare  their  gene  expression  in  LNCaP- Abl,  a  castration  resistant 
counterpart.  (Months  5-6) 

Verify  the  expression  of  PLZF  and  its  targets  in  LAPC4-PR/CR  and  VCaP/VCS2,  two  other  pairs 
of  hormone  naive  and  castration  resistant  cells.  (Months  7-8) 

b.  Integrated  analysis  of  PLZF  cistrome  and  PLZF  binding  regions  by  using  the  chromatin 
immunoprecipitation  sequencing  (ChIP-seq)  technique.  Perform  PLZF  ChIP-seq  in  the  presence 
and  absence  of  DHT  in  LNCaP.  (Months  9-10)  Bioinformatic  integration  of  PLZF  ChIP-seq  and 
microarray  data.  (Months  11-12) 

Task  2.  Investigate  the  biological  function  of  PLZF  in  CRPC  development.  (Months  13-18) 

a.  Determine  whether  FLZF-depletion  confers  prostate  cancer  cell  castration  resistant  growth  in  vitro. 
Perform  cell  proliferation  assay  in  reduced  medium  (C-FBS)  using  transient  RNAi  silencing  of 
PLZF  in  LNCaP  (Months  13-14) 

b.  Determine  whether  FLZF-depletion  promotes  CRPC  tumor  growth  in  vivo.  Generate  PLZF 


stable  shRNA  knockdown  cells.  (Months  15-18)  Conduct  mouse  xenograft  studies  using 
PLZF  stable  knocking  down  cells.  (19-24) 


What  was  accomplished  under  these  goals? 
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Figure  1:  Systematic  RNAi  suppressor  screen  in  androgen-sensitive  LNCaP  cells 

(a)  Left:  A  flowchart  of  screening  for  the  clones  showing  androgen-independent 
growth  using  a  pooled  genome-scale  shRNA  library.  Right:  Venn  diagram  showing 
overlapping  between  the  top  50  candidate  screen  hits,  PLZF  somatic  status  and  the  top 
60  DHT-induced  genes  (>2-fold).  (b)  The  colonies  were  stained  by  crystal  violet  (CV) 
and  photographed.  The  efficiency  and  efficacy  PLZF  shRNA  knockdown  was 
measured  by  Western  blot  (Insert).  CV  Staining  was  used  to  measure  the  cell  viability, 
(c)  RT-qPCR  and  Western  blotting  were  used  to  measure  PLZF  mRNA  and  protein 
expression  with  or  without  lOnM  of  DHT  and/or  lOpM  of  bicalutamide  treatment  in 
LNCaP  cells.  Values  were  the  mean  ±  SD  (n>3)  and  <  0.05. 


growth  suppressive  and  may  be  a  part  of  an  AR-repressed  oncogenic  program 
to  ADT. 


First,  in  support  of  a  role 
of  PLZF  in  the  CRPC  phenotype, 
we  used  an  androgen-dependent 
LNCaP  cell  line  as  a  model  to 
conduct  a  “Genome-Wide  RNAi 
Suppressor  Screen”,  aiming  to 
uncover  genes  whose  silencing 
were  crucial  for  androgen- 
independent  growth  (Botta  et  al., 
manuscript  in  preparation). 
Intriguingly,  among  the  top  50 
candidates  from  this  global  RNAi 
screen,  we  found  that  PLZF  was 
the  only  gene  that  was  also 
induced  by  androgens  (Fig.  la,  b, 
and  c),  suggesting  that  PLZF  is 
both  androgen-responsive  and 
essential  for  developing  resistance 


Next,  we  examined  the  growth  of  LNCaP  cells  using  2  additional  shRNAs  directed  at  PLZF.  These 
shRNAs  also  induced  androgen-independent  growth  in  LNCaP  cells;  whereas  non-specific  silenced  (shCtrl) 
cells  did  not  grow,  supporting  the  findings  of  the  screen  (Fig.  2a).  Interestingly,  while  androgen  depletion  alone 
(cultured  in  CSS)  lowers  PLZF  expression,  it  is  actually  the  further  reduction  of  PLZF  by  the  shRNAs  that 
promotes  growth. 
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Figure  2:  The  effect  of  PLZF  alternation 
in  prostate  cancer  cells  (a)  Growth  and  WB 
of  LNCaP  using  two  additional  PLZF 
shRNAs  (#1  and  2)  or  scrambled  control 
(shCtrl)  in  5%  charcoal-stripped  serum 
(CSS),  (b)  The  colonies  with  or  without 
PLZF  knockdown  were  cultured  in  (a) 
regular  or  (b)  CSS  medium  and  treated  with 
or  without  enzalutamide.  (d)  Growth  of 
LNCaP  stable  cell  lines  was  determined  by 
crystal  violet  assay  after  3  weeks.  WB  of 
LNCaP  cells  expressing  GFP  or  PLZF  ORF. 
Each  column  was  relative  to  the 
corresponding  the  first  columnand  shown  as 
mean  ±  SD  (n  >  3),  and  *p  <  0.05. 


To  block  potential  AR  action,  we  used  enzalutamide  to  examine  the  cell  growth.  Significantly,  PLZF- 
depleted  cells  initially  responded  to  enzalutamide  but  rapidly  developed  resistance,  regardless  of  culturing 


conditions,  indicating  that  LNCaP  cells  with  PLZF  inactivation  is  growth  promoting  even  in  the  absence  of 
androgen  and  further  when  blocked  by  enzalutamide  (Fig.  2b  and  c). 

As  a  complementary  approach,  we  also  addressed  whether  PLZF  has  tumor  suppressing  effects  on  cell 
viability.  Importantly,  overexpression  of  PLZF  resulted  in  a  profound  growth  inhibitory  effect  on  PLZF-depleted 
cells  as  early  as  at  the  first  week  (Fig.  2d).  It  should  be  noted  that  the  pictures  of  crystal  violet  (CV)  staining 
were  taken  after  3  weeks  of  CSS  incubation  which  allows  shPLZF  (#7816)/plx-GFP  colonies  to  accumulate. 
Taken  together,  our  in  vitro  experiments  functionally  validate  the  tumor  suppressive  role  of  PLZF  identified 
from  the  systematic  RNAi  suppressor  screen  and  show  that  PLZF  is  an  androgen-responsive  gene  involved  in 
growth  suppression. 

To  determine  the  requirement  of  PLZF 
loss  for  androgen-independent  growth  in  vivo ,  we 
conducted  a  mouse  xenograft  study  and  showed 
that  tumors  with  PLZF  depletion  had  a  substantial 
growth  advantage  as  compared  to  the  control 
(Ctrl)  arm  in  castrated  nude  mice  (Fig.  3a).  To 
broaden  our  in  vitro  and  in  vivo  findings,  we 
investigated  PLZF  gene  expression  in  different 
stages  of  human  prostate  tumors.  The  gene 
expression  of  PLZF  is  significantly  decreased  in 
mCRPC  compared  to  its  expression  in  primary 
prostate  tumors.  (Fig.  3b),  indicating  that  PLZF 
functions  as  a  tumor  suppressor  in  prostate  cancer 
and  its  gene  expression  is  decreased  in  mCRPC. 
Collectively,  our  in  vitro  and  in  vivo  together  with 
human  genomic  sequencing  data  strongly  support 
that  PLZF  functions  as  a  tumor  suppressor  in 
prostate  cancer  and  demonstrate  that  PLZF  inactivation  as  an  important  molecular  event  for  CRPC 
development. 

To  interrogate  the  biological  consequences  of  the  loss  of  PLZF  expression,  we  performed  a 
bioinformatic  analysis  of  PLZF  ChIP  seq  and  PLZF-regulated  gene  expression  profiling  to  define  the  PLZF 
transcriptional  program.  Since  AR  exerts  both  growth-promoting  and  growth-suppressing  functions  in 
maintaining  the  equilibrium  between  cell  differentiation  and  proliferation^],  we  hypothesized  that  ADT 
designed  to  disrupt  AR  signaling  may  stimulate  an  AR-repressed  oncogenic  program  through  PLZF  in  adaption 
to  ADT,  contributing  to  the  castration  resistant  phenotype  [9].  Interestingly,  we  uncovered  that  the  genes  whose 
expression  were  up-regulated  in  PLZF-depleted  cells  were  highly  associated  with  the  MAPK  pathway,  including 
5  PLZF  direct  targets,  RRAS,  MKNK2,  DDIT3,  JUND  and  JUN  (Fig.  4),  suggesting  that  PLZF  might  be 
regulating  MAPK  signaling  pathway.  Taken  together,  our  data  suggest  that  ADT  results  in  down-regulation  of 
PLZF,  which,  in  turn,  activates  MAKP  activity.  These  findings  are  consistent  with  the  findings  that  MAPK 
signaling  is  up-regulated  in  CRPC  patients  and  murine  models [10,  1 1].  In  summary,  we  propose  an  intrinsic 
resistance  mechanism  through  PLZF  down  regulation  or  loss,  wherein  an  AR-repressed  oncogenic  program 
allows  residual  prostate  tumor  cells  to  adapt  to  castrate  levels  of  androgens. 


Figure  3:  PLZF  functions  as  a  tumor  suppressor  in  vivo  (a) 

Tumor  formation  assays  of  castrated  male  nude  mice  injected  with 
shCtrl  and  PLZF  stable  silencing  22Rvl  cells.  IHC  of  PLZF  and  Ki- 
67  were  used  to  monitor  for  effects  and  efficacy  of  shPLZF  and  cell 
proliferation  in  22Rvl  xenografts.  The  photo  shows  each  22Rvl- 
shPLZF  tumor  relative  to  the  corresponding  control  on  the  opposite 
flank,  (b)  the  expression  of  PLZF  in  27  hormone-sensitive  prostate 
cancer  (HSPC)  and  29  mCRPC  tumors. 


What  opportunities  for  training  and  professional  development  has  the  project  provided? 

"Nothing  to  Report" 

How  were  the  results  disseminated  to  communities  of  interest? 

"Nothing  to  Report" 


What  do  you  plan  to  do  during  the  next  reporting  period  to  accomplish  the  goals? 

First,  I  will  validate  PLZF  recruitment  to  target  loci  and  regulation  on  their  gene  expression,  using  ChIP-  and 
RT-qPCR.  Next,  I  will  examine  PLZF-mediated  MAPK  pathway  activity  and  other  phosphoprotein  signaling. 
Second,  I  will  assess  the  effect  of  relevant  pharmacological  compounds  on  PLZF-ablated  prostate  cancer  cells. 
Initially,  I  will  concentrate  on  using  MAPK  inhibitors,  since  Dr.  Tom  Roberts’  laboratory  has  previously 
reported  that  the  MAPK  pathway  is  aberrantly  activated  in  the  murine  GEM  CRPC  model  in  vivo  [11]. 

IMPACT: 

Promyelocytic  Zinc  Finger  Protein  (PLZF)  is  genetically  lost  in  5-10%  of  tumors  from  metastatic  castration 
resistant  prostate  cancer  (mCRPC)  patients  who  previously  were  treated  with  androgen  deprivation  therapy 
(ADT)  and  developed  resistance.  In  addition,  PLZF  expression  can  be  down-regulated  by  ADT,  since  PLZF  is  a 
well-known  androgen  receptor  (AR)-stimulated  gene.  Our  data  demonstrated  that  prostate  cancer  cells  in  which 
PLZF  loss  is  mimicked  by  short  hairpin  RNA-mediated  knockdown  exhibited  a  CRPC  and  enzalutamide  (a 
second-generation  antiandrogen)  resistant  phenotype  in  vitro  and  in  vivo.  The  hope  is  that  this  approach  to 
define  that  patients  with  PLZF  somatic  deletions  represent  a  distinct  molecular  subtype  in  genetically 
heterogeneous  CRPC,  serving  as  a  potential  predicative  biomarker  of  drug  response  to  AR-directed  therapies,  as 
well  as  providing  a  basis  for  clinical  development  of  personalize  therapy.  As  such,  this  proposal  addresses  the 
Overarching  Challenge  of  “Developing  effective  treatments  and  addressing  mechanisms  of  resistance  for  men 
with  high-risk  or  metastatic  prostate  cancer”.  In  so  doing,  our  proposal  also  addresses  the  Focus  Areas  of 
“Genetics”,  “Mechanisms  of  Resistance”  and  “Therapy”. 


What  was  the  impact  on  the  development  of  the  principal  discipline(s)  of  the  project? 

Whole  exome  sequencing  of  metastatic  castration  resistant  prostate  tumors  (mCRPC)  from  the  Michigan  cohort 
and  the  East  Coast  Stand  Up  to  Cancer/Prostate  Cancer  Foundation  (SU2C/PCF)  study  recently  revealed  that 
5-10%  of  tumors  harbor  promyelocytic  leukemia  zinc  finger  (PLZF)  focal  homozygous  deletions.  This  cohort  of 
mCRPC  patients  who  were  previously  were  treated  with  androgen  deprivation  therapy  (ADT)  developed  castrate 
resistance.  Intriguingly,  from  a  “Genome-Wide  RNAi  Suppressor  Screen”,  aiming  to  uncover  genes  whose 
silencing  were  crucial  for  androgen-independent  growth  global  RNAi  screen,  we  found  that  PLZF  was  the  only 
gene  that  was  also  induced  by  androgens,  suggesting  that  PLZF  is  both  androgen-responsive  and  growth 
suppressive  and  may  be  a  part  of  an  AR-repressed  oncogenic  program  essential  for  developing  resistance  to  ADT. 


What  was  the  impact  on  other  disciplines? 

"Nothing  to  Report" 

What  was  the  impact  on  technology  transfer? 

"Nothing  to  Report" 

What  was  the  impact  on  society  beyond  science  and  technology? 

"Nothing  to  Report" 

CHANGES/PROBLEMS: 

"Nothing  to  Report" 

PRODUCTS: 

1.  Selected  poster  “Intensive  androgen  deprivation  therapy  compensatorily  activates  PLZF-repressed 
oncogenic  circuitry  that  reprograms  the  residual  prostate  cancer  “  for  the  AACR/PCF  Conference: 
Advances  in  prostate  cancer  research  (January  18-21,  2014,  San  Diego,  CA) 

2.  AACR-Alfac,  Incorporated  Scholar-in-Training  Award  by  AACR  (January,  2014) 

3.  Invited  short  talk  “Androgen  deprivation  therapy  activates  PLZF-repressed  oncogenic  circuitry  that 
reprograms  the  residual  prostate  cancer”  at  Seventh  Annual  Prostate  Cancer  Program  Retreat  (March  16- 
18,  2014,  Fort  Lauderdale,  FL) 

4.  Publication  in  a  peer-reviewed  journal: 

Hsieh  CL,  Teng  Fei,Yiwen  Chen,  Tiantian  Li,  Yanfei  Gao,  Xiaodong  Wang,  Tong  Sun,  Christopher  L. 
Sweeney,  Gwo-Shu  Mary  Lee,  X.  Shirley  Liu,  Myles  Brown  and  Philip  W.  Kantoff.  Enhancer  RNAs 
participate  the  androgen  receptor-driven  looping  that  selectively  enhances  gene  activation.  Proc  Natl 
Acad  Sci  USA.  2014  May  20;lll(20):7319-24. 


PARTICIPANTS  &  OTHER  COLLABORATING  ORGANIZATIONS 
What  individuals  have  worked  on  the  project? 
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The  androgen  receptor  (AR)  is  a  key  factor  that  regulates  the 
behavior  and  fate  of  prostate  cancer  cells.  The  AR-regulated  network 
is  activated  when  AR  binds  enhancer  elements  and  modulates 
specific  enhancer-promoter  looping.  Kallikrein-related  peptidase  3 
(KLK3),  which  codes  for  prostate-specific  antigen  (PSA),  is  a  well- 
known  AR-regulated  gene  and  its  upstream  enhancers  produce  bi¬ 
directional  enhancer  RNAs  (eRNAs),  termed  KLK3e.  Here,  we  dem¬ 
onstrate  that  KLK3e  facilitates  the  spatial  interaction  of  the  KLK3 
enhancer  and  the  KLK2  promoter  and  enhances  long-distance  KLK2 
transcriptional  activation.  KLK3e  carries  the  core  enhancer  element 
derived  from  the  androgen  response  element  III  (ARE  III),  which  is 
required  for  the  interaction  of  AR  and  Mediator  1  (Medl).  Further¬ 
more,  we  show  that  KLK3e  processes  RNA-dependent  enhancer  ac¬ 
tivity  depending  on  the  integrity  of  core  enhancer  elements.  The 
transcription  of  KLK3e  was  detectable  and  its  expression  is  signifi¬ 
cantly  correlated  with  KLK3  (R2  =  0.6213,  P  <  5  x  KT11)  and  KLK2 
(R2  =  0.5893,  P  <  5  x  1(T10)  in  human  prostate  tissues.  Interestingly, 
RNAi  silencing  of  KLK3e  resulted  in  a  modest  negative  effect  on  pros¬ 
tate  cancer  cell  proliferation.  Accordingly,  we  report  that  an  andro¬ 
gen-induced  eRNA  scaffolds  the  AR-associated  protein  complex  that 
modulates  chromosomal  architecture  and  selectively  enhances  AR-de- 
pendent  gene  expression. 

KLK3e/AR/Med1  complex  |  chromosomal  looping 

Androgens  and  the  androgen  receptor  (AR)  play  pivotal  roles 
not  only  in  the  development  of  the  normal  prostate  gland 
but  also  in  the  growth  and  progression  of  prostate  cancer  (PCa) 
(1-4).  Androgens  exert  their  biological  and  physiological  effects 
through  activating  AR  transcriptional  activity.  AR,  a  transcrip¬ 
tion  factor,  recruits  ligand-dependent  transcriptional  machinery 
that  activates  a  tissue-specific  transcriptional  program  involved 
in  development  and  differentiation  (4-6).  In  PCa  cells,  AR 
occupies  numerous  gene  loci  to  activate  transcripts  mediating 
critical  cellular  activities,  such  as  energy  and  metabolism  (5).  AR 
activation  and  reactivation  are  essential  for  PCa  development 
and  tumor  recurrence  (7). 

One  unanticipated  observation  that  emerged  from  multiple 
global  studies  of  RNA  polymerase  II  (Pol  II)  occupancy  was  that 
Pol  II  binds  to  a  large  number  of  intergenic  AR-bound 
enhancers  marked  by  histone  H3  lysine  4  monomethylation 
(H3K4mel)  and  lysine  27  acetylation  (H3K27ac)  and  produces 
enhancer-derived  long  noncoding  RNAs  (eRNAs)  (8-11).  Al¬ 
though  eRNAs  are  bidirectional  and  nonpolyadenylated,  their 
induction  correlates  with  the  induction  of  adjacent  exon-coding 
genes  (9).  The  activity  of  eRNAs  was  shown  to  augment  the 
expression  of  not  only  neighboring  but  also  intrachromosomally 
distant  genes  to  convey  enhancer-dependent  transcription  (12, 
13).  Recently,  the  5'-GRO-seq  data  have  revealed  that  sequence- 
specific  transcription  factors  and  eRNA  transcripts  are  required 
for  transfer  of  full  enhancer  activity  to  a  target  promoter  (14). 


Moreover,  eRNAs  have  been  shown  to  configure  a  chromatin 
state  that  facilitate  transcription  factors  binding  on  the  target 
promoters  at  defined  genomic  loci  (15). 

Mechanistically,  eRNAs  have  been  shown  to  associate  with  the 
cohesin  complex  at  estrogen  receptor  binding  sites  and  contrib¬ 
ute  to  gene  activation  by  stabilizing  estrogen  receptor-a  (ER- 
a)/eRNA-induced  enhancer-promoter  looping  (16,  17).  In  ad¬ 
dition,  ncRN A- activating  (ncRNA-a),  a  class  of  enhancer-like 
RNAs,  has  been  shown  to  interact  with  the  Mediator  complex 
that  enhances  transcription  in  a  c/s-mediated  mechanism  (18), 
highlighting  an  RNA-dependent  functional  effect.  Therefore,  we 
hypothesized  that  eRNAs  facilitates  spatial  communication  of 
enhancers  and  promoters,  thereby  enhancing  transcriptional 
activation  events.  Here,  we  report  that  the  Kallikrein-related 
peptidase  3  (KLK3)  enhancer,  one  of  the  strongest  AR-bound 
enhancers  in  prostate  cancer  cells,  produces  KLK3  eRNA 
(KLK3e),  which  impacts  androgen-induced  gene  activation.  This 
work  reveals  that  KLK3e  selectively  enhances  AR-regulated 
gene  expression  through  its  collaboration  with  the  AR-de- 
pendent  looping  complex.  Thus,  our  findings  provide  evidence 
for  a  positive  regulatory  loop  in  which  the  AR-regulated  network 
is  modulated  by  an  intermediate  eRNA. 

Results 

The  KLK3  Enhancer  Produces  eRNA  (KLK3e).  We  compiled  recently 
established  genomic  datasets  (19)  and  showed  that  the  KLK3 
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We  report  that  enhancer  RNAs  (eRNAs),  a  class  of  long  non¬ 
coding  RNAs,  participate  in  the  androgen  receptor  (AR)- 
dependent  looping  complex  that  enhances  spatial  communica¬ 
tion  of  distal  enhancers  and  target  promoters,  leading  to  tran¬ 
scriptional  activation  events.  Furthermore,  our  data  show  that 
KLK3  eRNA  (KLK3e)  selectively  enhances  the  gene  expression  of 
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durability  of  androgen  deprivation  therapy  (ADT)  and  prediction 
of  the  efficacy  of  ADT  by  measuring  the  enhancer-derived  ac¬ 
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enhancer,  also  defined  as  the  androgen  response  element  III 
(ARE  III)  (20),  was  marked  by  AR,  H3K27ac,  and  H3K4mel, 
implying  that  this  region  is  transcriptionally  active  (Fig.  L4).  In 
fact,  the  transcription  of  KLK3  eRNA  (KLK3e)  was  demon¬ 
strated  by  GRO-seq  data  (10)  and  confirmed  by  our  reverse 
transcription-quantitative  PCR  (RT-qPCR)  analyses  in  andro¬ 
gen-dependent  LNCaP  and  VCaP,  and  androgen-independent 
LNCaP-abl  cells  (21)  (Fig.  IB  and  Fig.  SI  A  and  B).  The  ex¬ 
pression  of  KLK3e  was  induced  within  the  first  hour  of  dihy¬ 
drotestosterone  (DHT)  exposure  and  its  expression  level  gradually 
increased  in  a  time-dependent  manner,  resembling  the  expression 
of  mature  KLK3  in  both  FNCaP  and  VCaP  cells;  in  contrast,  al¬ 
though  transcription  at  the  KLK3  locus  remains  active  in  andro¬ 
gen-independent  FNCaP-abl,  androgens  (DHT)  did  not  stimulate 
their  mRNA  expression  to  a  comparable  magnitude  as  those  in 
parental  FNCaP  cells  (Fig.  SIB).  Next,  we  investigated  the  impact 
of  AR  antagonist,  bicalutamide,  on  the  expression  of  KLK3e.  The 
expression  of  KLK3e  and  KLK3  was  induced  upon  DHT  stimu¬ 
lation  and  blocked  by  bicalutamide  (Fig.  SIC),  revealing  that  li¬ 
gand-dependent  eRNA  synthesis  takes  place  in  association  with 
target  mRNA  transcription.  Taken  together,  our  findings  are 
consistent  with  the  previous  report  that  there  is  a  high  correlation 
of  activity-dependent  induction  between  eRNAs  and  adjacent 
coding  mRNAs  (9). 
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Fig.  1.  The  KLK3  enhancer  produces  IncRNAs.  (A)  The  chromatin  status  in 
the  vicinity  of  KLK3.  ChIP-seq  datasets  of  AR,  H3K27ac,  and  H3K4me1  (19). 
(B)  RT-qPCR  analyses  of  expression  of  antisense  KLK3  eRNA  (red  bar),  sense 
KLK3  eRNA  (green  bar),  and  mature  KLK3  (blue  bar)  upon  DHT  (10  nM) 
treatment  in  a  time-dependent  manner  in  LNCaP  cells.  (C,  Upper)  A  diagram 
of  functional  ARE  positions  that  are  relative  to  the  KLK3  transcription 
starting  site  (TSS)  adapted  from  ref.  20.  (Lower)  RT-qPCR  analyses  of  the 
expression  of  segments  of  KLK3e  with  DHT  (10  nM)  treatment  in  LNCaP  cells. 
The  regions  of  LNA  probes  against  the  antisense  and  sense  KLK3e  were  as 
indicated.  Primer  sets  for  qPCR  were  designed  to  specifically  target  either 
antisense  (primers  AS1-AS4)  or  sense  strand  (primer  S1-S7)  of  KLK3e  as 
shown.  The  right  electrophoresis  images  show  the  PCR  products  of  antisense 
(2,230  bp)  and  sense  (2,170  bp)  KLK3e.  All  data  analyzed  by  RT-qPCR  were 
normalized  to  the  expression  level  in  vehicle  (-).  Data  are  shown  as  mean  ± 
SD  (n  >  3). 


Although  multiple  studies  have  identified  different  eRNA 
species  from  mammalian  enhancers,  the  precise  mapping  of 
eRNAs  has  not  been  described.  Because  eRNAs  produced  from 
these  intergenic  enhancers  are  low  at  steady  state  (9)  and  esti¬ 
mated  at  fewer  than  100  copies  per  cell  (17),  we  enriched  cellular 
RNAs  by  depleting  rRNAs  before  reverse  transcription.  PCR 
primer  sets  against  the  antisense  (AS1-AS4)  and  sense  (S1-S7) 
of  KLK3e  were  designed  based  on  the  GRO-seq  and  genomic 
datasets  (10,  19).  As  shown  in  Fig.  1C,  the  transcription  of  bi¬ 
directional  KLK3e  was  captured;  whereas,  no  robust  DHT- 
induced  transcription  was  detectable  beyond  AS4  and  S7.  We 
performed  primer-extension  PCR  assay  to  connect  shorter  frag¬ 
ments  (Fig.  SID)  and  further  showed  PCR  products  for  the  an¬ 
tisense  (2,230  bp)  and  sense  (2,170  bp)  of  KLK3e  (Fig.  1C  and 
cDNA  sequences  available  in  SI  Experimental  Procedures).  To  dis¬ 
tinguish  whether  the  eRNA  at  KLK3  enhancer  are  polyadenylated 
or  not,  we  performed  reverse  transcription  using  random  hexamer 
primers  and  3 '-RACE  (poly-T)  adapter  from  the  same  RNA 
materials.  As  shown  in  Fig.  SIT,  the  expression  of  AS4  and  S7  was 
substantially  induced  by  DHT  treatment,  suggesting  that  the  poly-A 
tail  was  added  to  the  3 '-end  of  antisense  and  sense  KLK3e.  In 
contrast,  SI  showed  less  DHT-stimulated  expression  in  the  poly-T 
priming  compared  with  random  hexamer  primed  cDNA.  These 
data  suggest  that  this  locus  produces  a  mixture  of  RNA  species; 
some  of  them  are  polyadenylated,  whereas  others  may  be  actively 
undergoing  synthesis  of  polyadenylation,  but  we  cannot  totally  ex¬ 
clude  the  existence  of  those  species  without  poly-A  tails.  To  detect 
the  size  of  the  mature  eRNA  at  this  loci,  we  customized  double 
DIG-labeled  probes  using  locked  nucleic  acids  (FNA)  (Exiqon). 
Northern  analysis  revealed  that  a  major  form  (~2.5  kb)  and  a  minor 
form  (~1.7  kb)  of  sense  KLK3e  were  produced  upon  DHT  in¬ 
duction,  whereas  no  signal  was  detected  at  the  antisense  strand  (Fig. 
SIT),  suggesting  that  the  expression  of  antisense  strand  may  be 
lower  than  the  levels  that  can  be  detected  by  the  FNA  probe  (Fig. 
SIT).  Collectively,  this  suggests  that  the  sense  strand  (S1-S7) 
KLK3e,  a  ~2.2-kb  transcript,  is  polyadenylated  and  yields  a  ~2.5-kb 
transcript  shown  in  Northern  analysis.  Although  the  antisense 
KLK3e  can  be  polyadenylated  and  detected  by  qPCR,  its  expression 
may  be  substantially  lower  than  the  sense  strand  and  cannot  be 
detected  by  Northern  analysis  under  the  same  experimental  settings, 
indicating  that  the  sense  strand  KLK3e  is  the  dominant  species. 

The  Regulatory  Function  of  KLK3e  on  the  Kallikrein  Locus.  To  in¬ 
vestigate  the  potential  function  of  eRNAs  in  transcriptional  reg¬ 
ulation,  short  interfering  RNAs  (siRNAs)  targeting  KFK3e 
(siKLK3e)  were  designed  on  the  basis  of  KLK3e  mapping. 
Results  obtained  from  distinct  sets  of  siRNAs  against  the  sense 
strand  of  KLK3e  showed  that  siKLK3e  specifically  targeted 
sense  KFK3e  and  significantly  impeded  DHT-induced  KLK3 
expression  in  both  FNCaP  and  VCaP  cells  but  had  a  minimal 
impact  on  FNCaP-abl  cells  (Fig.  2 B  and  Fig.  S2  A-E),  consistent 
with  the  notion  that  ligand-activated  eRNAs  are  important  for 
nearby  gene  activation  (13,  17).  Considering  that  KLK15,  2,  and 
4  were  adjacent  to  KFK3  and  androgen-responsive  (Fig.  24),  we 
examined  whether  KFK3e  impacted  transcription  of  these 
KLKs.  Interestingly,  knocking  down  of  sense  KLK3e  significantly 
inhibited  the  expression  of  KLK2,  whereas  other  KLKs  remained 
unchanged  (Fig.  2 B).  We  also  investigated  the  effect  of  antisense 
KLK3e  on  the  regulation  of  target  gene  expression.  RNAi  si¬ 
lencing  of  the  antisense  of  KLK3e  showed  a  marginal  inhibitory 
effect  on  KLK3  and  KLK2  gene  expression  (Fig.  S2 T).  Despite 
the  presence  of  bidirectional  KLK3e  transcription,  data  from 
siRNA  to  KLK3e  and  Northern  analysis  showed  that  the  sense 
strand  of  KLK3e  is  the  dominant  transcript  with  functional  effects  on 
the  expression  of  KLK3  and  distant  KLK2.  Thus,  we  chose  the  sense 
KLK3e  and  continued  to  study  its  RNA-dependent  enhancer  activity. 

Because  the  induction  of  eRNA  is  strongly  correlated  with  the 
enhancer  looping  to  target  gene  promoters  (13,  16,  17),  it  is 
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Fig.  2.  The  regulatory  function  of  KLK3e  on  the  KLK  locus.  (A)  A  heat  map 
of  the  KLK  transcripts  upon  4  and  16  h  of  DHT  exposure  in  LNCaP  cells  (21). 
Genes  are  listed  as  their  order  in  the  kallikrein  loci  from  centromere  (Top)  to 
telomere  (Bottom).  (B,  Upper)  Schematic  drawing  indicates  the  relative  po¬ 
sition  of  AR-regulated  KLK  genes  and  KLK3e  on  the  chromosome  1 9.  (Lower) 
LNCaP  cells  were  treated  with  vehicle  (-)  or  DHT  (10  nM)  and  50  nM  of 
nonsilencing  control  (siCtrl)  or  KLK3e  siRNA  (siKLK3e);  lowercase  "e"  deno¬ 
ted  as  eRNA.  Efficacy  and  effects  of  KLK3e  depletion  on  KLK  genes  were 
assessed  by  RT-qPCR.  (C,  Upper)  Schematic  graph  shows  the  KLK3/2  gene 
loci.  (Lower)  LNCaP  cells  were  treated  with  vehicle  or  DHT  for  4  h.  KLK2 
mRNA  expression  was  assessed  by  RT-qPCR  (Inset).  The  relative  cross-linking 
frequency  between  the  anchor  region  (ARE  III)  and  distal  fragments  (shaded 
bars)  was  determined  by  qPCR  and  normalized  to  the  control  region  (frag¬ 
ment  A).  (D)  The  relative  cross-linking  frequency  between  the  anchor  region 
D  and  E  region  was  determined  by  qPCR  and  normalized  to  the  A  region.  All 
data  analyzed  by  RT-qPCR  or  qPCR  were  normalized  to  the  expression  level 
in  vehicle  and  are  shown  as  mean  ±  SD  (n  >  3),  and  *P  <  0.05. 


possible  that  the  KLK3  enhancer  element  is  able  to  interact  with 
the  KLK2  promoter.  To  address  this,  we  used  the  chromosome 
conformation  capture  (3C)  assay  to  study  the  spatial  interaction 
of  the  KLK3/2  loci.  It  should  be  noted  that  the  KLK2  promoter 
contains  an  ARE  I  that  shares  over  80%  sequence  homology 
with  the  ARE  I  in  the  KLK3  promoter  (20).  Indeed,  DHT 
induces  KLK2  transcription  and  the  interaction  between  the 
anchor  and  fragment  D  and  E  by  twofold  to  threefold  (Fig.  2 C). 
Amplified  PCR  of  anchor-fragment  D  (105  bp)  and  anchor- 
fragment  E  (115  bp)  was  confirmed  by  DNA  sequencing  (Fig.  S2 
G  and  H),  suggesting  that  the  KLK3  enhancer  interacts  with  the 
KLK2  promoter  to  mediate  distal  KLK2  gene  activation. 

To  study  the  mechanism  of  KLK3e-mediated  transcriptional 
events,  we  looked  for  KLK3e  interacting  partners  by  asking 
whether  KLK3e  interacts  with  AR,  because  KLK3e  carries  the 
core  enhancer  element.  We  conducted  AR-directed  RNA  im- 
munoprecipitation  (RIP)  coupled  with  qPCR  using  primer  sets 
targeting  the  core  enhancer  element  at  5'-end  (SI),  middle 
segment  (S4),  and  3'-end  (S7)  of  the  sense  KLK3e.  Our  results 
show  a  specific  DHT-induced  interaction  between  AR  and  SI 


(Fig.  S27),  suggesting  that  AR  relies  on  the  core  enhancer  ele¬ 
ment  of  KLK3e.  Because  Medl  has  been  shown  to  directly  in¬ 
teract  with  AR  and  is  involved  in  the  chromosomal  interaction  of 
enhancer-promoter  looping  (22,  23),  we  next  asked  whether 
Medl  participates  in  the  observed  eRNA  functions.  RIP  analysis 
revealed  a  strong  androgen-stimulated  interaction  of  SI  and 
Medl  (Fig.  S27),  comparable  to  the  interaction  of  KLK3e  and 
AR.  Conversely,  the  antisense  KLK3e  (AS1)  exhibited  modest 
interaction  with  AR  and  Medl  (Fig.  S2 K),  implying  that  the 
enhancer  function  of  these  two  KLK3e  species  is  determined  by 
the  levels  of  interaction  with  the  AR-Medl  complex. 

To  investigate  the  potential  role  of  KFK3e  in  the  KLK3/2  loci, 
we  next  used  3C-qPCR  to  assess  the  interaction  of  the  anchor- 
fragment  D/E  regions  in  KLK3e-depleted  cells.  Silencing  of 
KLK3e  impaired  KFK2  gene  expression  and  the  interaction 
of  the  KLK2  promoter  and  the  KLK3  enhancer  (Fig.  2D).  Sup¬ 
pression  of  Medl  also  substantially  reduced  the  interaction  of  the 
KLK3/2  loci  (Fig.  2D  and  Fig.  S2L),  supporting  a  pivotal  role  of 
Medl  in  mediating  this  long-range  chromatin  looping  (18,  22). 
Hence,  our  results  demonstrate  that  both  KLK3e  and  Medl  are 
involved  in  maintaining  the  spatial  interaction  between  an  en¬ 
hancer  and  a  target  promoter. 

The  KLK3e/AR/Med1  Ribonudeoprotein  Complex  Transcriptionally  Regulates 
Target  Promoters.  To  ensure  that  the  endogenous  KLK3e  plays  a 
functional  role  on  its  target  promoters,  we  performed  chromatin 
immunoprecipitation  (ChIP)  using  antibodies  against  AR  and  ac¬ 
tivated  RNA  polymerase  II  (Pol  II  S5p)  in  KLK3e-depleted  cells. 
These  analyses  demonstrated  that  depletion  of  KLK3e  inhibited 
AR  binding  and  Pol  II  activation  at  the  KLK2  promoter,  consistent 
with  the  reduction  of  gene  expression  of  KLK2  (Fig.  3 A).  Next,  we 
examined  the  regulatory  effect  of  Medl  on  the  target  promoters. 
As  expected,  knockdown  of  Medl  significantly  reduced  AR  oc¬ 
cupancy  and  Pol  II  activation  at  the  KLK2  promoter  (Fig.  3 B). 
Whether  KLK3e  collaborates  with  Medl  and  functionally  affects 
transcription  of  the  KLK2  promoter  is  unclear.  Depletion  KLK3e 
resulted  in  decrease  of  Medl  binding  to  the  promoters  of  KLK2 
(Fig.  3 C).  In  addition,  the  effect  of  siKLK3e  and  siMedl  on  the 
gene  expression  of  KLK2  and  Medl  was  shown  in  Fig.  S3A.  Al¬ 
together,  our  data  suggest  that  KLK3e  cooperates  with  Medl  to 
enhance  KLK2  gene  expression.  It  should  be  noted  that  the  actions 
of  AR,  Medl,  and  Pol  II  at  the  KLK3  promoter  were  comparable 
to  those  at  the  KFK2  promoter  (Fig.  S3  B-D ),  implying  that  the 
regulatory  function  of  KLK3e  at  the  KLK3/2  promoters  may  act 
through  the  same  AR-dependent  ribonudeoprotein  complex. 

Interestingly,  silencing  of  KLK3e  did  not  affect  the  action  of 
AR  and  Pol  II  on  the  KLK3  enhancer  (Fig.  S3 E),  suggesting  that 
the  biogenesis  of  KLK3e  takes  place  after  ligand-activated  AR 
was  recruited  to  the  enhancer.  Although  depletion  of  Medl 
blocked  AR  binding  and  Pol  II  activation  on  the  KFK3  en¬ 
hancer,  the  expression  of  KFK3e  was  not  significantly  affected, 
implying  that  residual  Medl -associated  module  may  be  sufficient 
for  KLK3e  production  (Fig.  S3F).  Collectively,  our  data  show  the 
mechanism  by  which  KLK3e  forms  a  functional  complex  with 
AR  and  Medl  to  promote  transcription,  supporting  the  findings 
that  transfer  of  full  enhancer  activity  to  a  target  promoter 
depends  on  enhancer-bound  transcription  factors  and  enhancer- 
derived  RNA  transcripts  (14). 

The  Core  Enhancer  Element  Renders  KLK3e  RNA-Dependent  Enhancer 
Activity.  To  validate  the  enhancer  function  of  KLK3e,  we  used 
a  luciferase  reporter  assay  to  quantitatively  measure  AR  tran¬ 
scriptional  activity.  Although  inclusion  of  KLK3e  (S1-S7)  or  flip¬ 
ped  (S7-S1)  sequences  markedly  increased  KLK2  promoter 
activity  by  fourfold  to  fivefold,  this  effect  was  abolished  when  ARE 
III  sequences  were  deleted  (S2-S7)  (Fig.  44),  confirming  the  im¬ 
portance  of  the  enhancer  integrity  and  their  orientational  inde¬ 
pendency  (17,  24). 
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Fig.  3.  The  KLK3e/AR/Med1  ribonucleoprotein  complex  transcriptionally 
regulates  target  promoters.  Chromatin  immunoprecipitation  (ChIP)  of  AR  or 
Medl  or  activated  Pol  II  (S5p)  was  performed  in  KLK3e-  or  Med  1 -depleted 
LNCaP  cells  with  or  without  DHT  (10  nM)  for  2  or  8  h.  (A-C)  qPCR  was 
conducted  to  measure  the  action  of  AR,  Pol  II  (S5p),  and  Medl  at  the  KLK2 
promoter.  Data  are  shown  as  mean  ±  SD  (n  >  3)  and  *P  <  0.05. 


The  expression  of  the  KLK3e  from  the  luciferase  reporter  was 
validated  by  RT-PCR  (Fig.  S44),  indicating  that  AREs  have  the 
capacity  to  produce  eRNAs.  We  next  asked  whether  the  KLK3e 
transcript  alone  is  adequate  for  this  transcriptional  enhance¬ 
ment.  For  this  purpose,  we  performed  a  cotransfection  study 
based  on  reporters  together  with  CMV-driven  sense  KLK3e 
constructs.  As  shown,  the  full-length  (S1-S7)  and  flipped  (S7-S1) 
KLK3e  increased  luciferase  expression  by  2-  -  2.5-fold,  whereas 
ARE  deletion  (S2-S7)  failed  to  enhance  promoter  activity  (Fig. 
4 B).  Furthermore,  antisense  KLK3e  (AS1-AS4)  was  unable  to 
induce  luciferase  expression  as  sense  KLK3e  (S1-S7)  did  (Fig.  S4 
B  and  C),  again  arguing  that  the  sense  KLK3e  is  the  functional 
long  noncoding  RNA  (IncRNA)  rather  than  the  antisense  strand. 
To  examine  whether  AR  is  required  for  the  action  of  KLK3e,  we 
conducted  experiments  in  AR-negative  COS-7  cells  and  confirmed 
that  the  activity  of  KLK3e  depends  on  the  availability  of  ligand- 
activated  AR  (Fig.  S4  D  and  E).  Altogether,  these  data  show  that 
the  core  enhancer  elements  are  essential  for  the  enhancer  activity 
of  KLK3e  and  that  ectopic  KLK3e  induces  the  promoter  activity 
on  a  luciferase  reporter,  implying  a  transregulatory  activation  me¬ 
diated  by  KLK3e. 

As  a  complementary  approach,  we  also  addressed  whether  the 
observed  luciferase  expression  is  mediated  through  endogenous 
KLK3e.  RNAi  silencing  of  KLK3e  significantly  reduced  lucif¬ 
erase  activity  in  the  reporters  with  or  without  KLK3e  insertion 
(Fig.  4C,  Left),  showing  that  KLK3e  is  required  for  the  enhance¬ 
ment  of  DHT-induced  activation.  When  KLK3e  was  ectopically 
expressed,  it  restored  the  luciferase  expression  inhibited  by  KLK3e 
suppression  (Fig.  S4F).  Most  importantly,  the  luciferase  activity  in 
the  reporter  with  KLK3e  insertion  was  reduced  by  -50%  in  Medl- 
depleted  cells  (Fig.  4 C,  Right),  suggesting  that  Medl  facilitates 
KLK3e-dependent  transcription.  Notably,  ncRNA-a  has  been 
shown  to  mediate  long-range  transcriptional  activation  through  its 
association  with  the  Medl  complex  (18).  In  summary,  our  results 


strongly  suggest  that  KLK3e  forms  a  functional  complex  with  AR 
and  Medl  that  facilitates  the  association  of  AR-bound  enhancers 
with  promoters  of  target  genes,  resulting  in  transcriptional  activation. 

KLK3  eRNA  Selectively  Enhances  AR-Regulated  Gene  Expression.  To 

explore  the  potential  effect  of  KLK3e  beyond  the  KLK  locus,  we 
measured  the  expression  of  canonical  androgen-regulated  genes 
in  KLK3e-depleted  cells.  Suppression  of  KLK3e  decreased  the 
expression  of  NKX3.1,  FKBP5,  and  PLZF,  whereas  the  expres¬ 
sion  of  TMPRSS2,  PDE9A,  and  GUCY1A3  was  not  affected 
(Fig.  5A).  Additionally,  we  observed  a  comparable  KLK3e- 
mediated  inhibition  of  NXK3.1  gene  expression  in  VCaP  cells 
(Fig.  S2D).  Moreover,  we  showed  the  expression  eRNAs  at  up¬ 
stream  of  NKX3.1,  FKBP5,  PLZF  was  induced  by  DHT  but  not 
affected  by  KLK3e  siRNA  knockdown  (Fig.  S5A),  suggesting  that 
the  KLK3e  siRNA  does  not  interfere  with  the  transcription  of 
eRNAs  found  at  the  enhancers  of  NKX3.1,  FKBP5,  and  PLZF. 
Although  eRNA-mediated  trans  effects  are  likely  to  be  relatively 
infrequent  or  quantitatively  small  (17),  our  data  show  that  KLK3e 
selectively  regulates  the  expression  of  AR-regulated  genes  that  do 
not  reside  on  the  same  chromosome,  implying  a  transregulatory 
function  of  KLK3e  on  AR-dependent  gene  expression. 

We  further  assessed  the  regulatory  role  of  KLK3e  and  Medl 
on  the  promoter  of  NKX3.1  and  TMPRSS2  to  mechanistically 
validate  our  above  findings  in  the  KLK  loci.  Although  depletion 
of  KLK3e  significantly  reduced  the  enrichment  of  AR  and  Pol  II 
S5p  at  the  NKX3.1  promoter  (Fig.  S 5B,  Upper),  the  action  of 
these  factors  at  the  TMPRSS2  promoter  was  unaffected  (Fig. 
S 5B,  Lower).  In  contrast,  Medl  knockdown  considerably  blocked 
both  AR  occupancy  and  Pol  II  activity  at  both  the  NKX3.1  and 
TMPRSS2  promoters  (Fig.  S5C).  These  results  demonstrate  that 
KLK3e  selectively  enhances  gene  expression,  whereas  Medl  has 
a  more  general  effect  on  the  regulation  of  AR-dependent  tran¬ 
scription.  Most  importantly,  KLK3e  depletion  reduced  Medl 
recruitment  to  the  promoter  of  NKX3.1,  not  TMPRSS2  (Fig.  S5 
D  and  E),  revealing  that  this  eRNA/AR/Medl  complex  may 
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Fig.  4.  The  core  enhancer  element  renders  KLK3e  RNA-dependent  en¬ 
hancer  activity.  (A,  Upper )  Schematic  representation  shows  the  insertion  in 
KLK2  promoter  driven  luciferase  reporter.  (Lower)  Reporter  assay  analyses 
were  performed  with  or  without  insertion  of  full-length  (S1-S7),  flipped  (57- 
SI),  and  ARE  deleted  (S2-S7)  of  sense  KLK3e.  ( B )  Reporter  analysis  of 
cotransfecting  the  reporter  with  independent  overexpressing  vector  as  in¬ 
dicated.  (C)  The  effects  of  siKLK3e  (Left)  or  si  Medl  (Right)  on  the  KLK2 
promoter-driven  Luciferase  activity.  Luciferase  activity  was  measured  in  the 
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Fig.  5.  KLK3e  selectively  enhances  AR-regulated  gene 
expression.  (A)  LNCaP  cells  were  treated  with  vehicle  or  DHT 
(10  nM)  and/or  siCtrl  or  siKLK3e.  The  effect  of  siKLK3e  on 
AR-regulated  gene  expression  was  determined  by  RT-qPCR. 
(B)  Correlation  between  KLK3e  (SI)  and  target  gene  ex¬ 
pression  in  vivo.  All  data  were  normalized  for  GAPDH  by 
the  differences  in  threshold  cycles  (ACT)  from  the  normal 
human  prostate  gland  (n  =  4),  primary  (n  =  11)  and  meta¬ 
static  (n  =  32)  prostate  tumors  (25).  (C)  LNCaP  cells  were 
cultured  in  5%  charcoal/dextran-treated  FBS  (C-FBS)  me¬ 
dium,  followed  by  siRNA  knockdown  of  scrambled  control 
(siCtrl)  or  KLK3e  (siKLK3e)  with  ethanol  (veh)  or  DHT  (10 
nM)  treatment.  Cell  growth  was  measured  by  WST1  assay. 
Data  represent  mean  ±  SD  (n  >  3),  and  *P  <  0.05.  (D)  A 
proposed  model  for  functional  activities  of  the  eRNA/AR/ 
Medl  complex  in  DHT-induced  c/s  (solid  line)  and  potential 
trans  (dish  line)  gene  activation  events. 


possess  a  degree  of  specificity  in  the  AR-dependent  transcrip¬ 
tional  machinery.  Furthermore,  we  evaluated  the  gene  expres¬ 
sion  from  the  normal  human  prostate  gland  (n  =  4),  primary  (n  = 
11)  and  metastatic  (n  =  32)  human  prostate  tumors  (25),  and 
showed  that  the  correlation  ( R 2)  of  KLK3e  (SI)  with  KLK3, 
KLK2,  NKX3.1,  and  TMPRSS2  was  0.6213  (P  <  5  x  KT11), 
0.5893  (P  <  5  x  1(T10),  0.4624  (P  <  2  x  KT7),  and  0.3787  (P  <  5  x 
KT6),  respectively.  In  addition,  the  expression  the  KLK3  and 
KLK2  was  the  most  correlated  (R2  =  0.8361,  P  <  3  x  10-16)  (Fig. 
5B  and  Fig.  S5F).  Collectively,  these  in  vivo  data  indicate  that  the 
transcription  of  KLK3e  was  detectable  and  its  presence  correlates 
with  KLK3/2  gene  expression.  We  noted  less  of  a  correlation 
between  KLK3e  expression  and  those  genes  in  trans ,  and  the 
correlation  was  stronger  with  NKX3.1  than  TMPRSS2,  consistent 
with  our  in  vitro  data,  although  these  data  do  not  prove  a  direct 
effect  in  vivo  of  KLK3e  on  the  transcription  of  these  genes.  Next, 
we  investigated  the  biological  effect  of  KLK3e  siRNA  on  LNCaP 
cells  and  showed  that  KLK3e  depletion  significantly  compromised 
cell  growth  in  the  presence  and  absence  of  exogenous  androgens 
(Fig.  5C).  Although  this  negative  effect  is  modest,  we  postulate 
that  suppression  of  multiple  eRNAs  produced  from  distinct  en¬ 
hancer  elements  may  synergetically  block  the  AR  transcriptional 
program.  Thus,  substantial  efforts  are  needed  to  fully  characterize 
the  transcriptional  and  biological  functions  of  other  significant 
androgen-induced  eRNAs  in  the  future. 

Discussion 

We  demonstrate  that  bidirectional  KLK3e  was  produced  from 
one  of  the  strongest  AR-bound  enhancers  (ARE  III)  and  its 
expression  pattern  resembles  KLK3  in  prostate  cancer  cell  lines 
(Fig.  1  and  Fig.  SI),  in  agreement  with  the  previous  report  that 
a  high  correlation  of  activity-dependent  induction  exists  between 
eRNAs  and  adjacent  coding  mRNAs  (9).  Suppression  of  KLK3e 
exclusively  reduces  gene  expression  of  KLK3  and  KLK2,  not  other 
KLK  family  members  and  that  KLK3e  facilitates  the  assembly 
of  the  transcriptional  apparatus  at  the  promoters  during  gene  ac¬ 
tivation  events.  This  classifies  KLK3e  as  a  family  of  functional 
IncRNA  involved  in  enhancing  long-distance  gene  transcription. 


Whether  eRNAs  function  as  molecular  bridges  that  mediate 
spatial  interactions  of  distal  enhancers  and  target  promoters,  or 
directly  configure  a  chromatin  state  that  facilitate  transcription 
factor  binding  at  the  target  promoters  is  unclear.  Recently, 
chromatin  isolation  by  RNA  purification  coupled  with  deep  se¬ 
quencing  was  deployed  to  interrogate  the  nature  of  RNA  and 
DNA  interaction  in  a  genome -wide  scale  (26).  Rosenfeld  and 
colleagues  (17)  used  this  methodology  to  identify  potential  sites 
where  estrogen-induced  FOXC1  eRNA  localizes  in  the  genome. 
Although  FOXC1  eRNA  is  specifically  enriched  in  a  small 
fraction  of  genomic  loci,  there  were  no  known  estrogen-regu¬ 
lated  genes  found  adjacent  to  these  RNA/DNA  binding  sites 
(17),  implying  that  the  functional  effect  of  eRNAs  on  target 
genes  may  not  completely  rely  on  nucleotide  base  pairing.  In  this 
report,  we  demonstrate  that  KLK3e  mediates  chromosomal 
looping  of  KLK3/2  loci  and  confirmed  in  a  heterologous  reporter 
assay  that  KLK3e  cooperates  with  AR  and  Medl  to  enhance 
KLK2  promoter  activity  through  the  core  enhancer  element 
(ARE  III).  The  state  of  KLK3e  and  KLK3/2  was  further  vali¬ 
dated  in  human  prostate  tissues  (n  =  47)  and  showed  that  the 
expression  of  KLK3e  significantly  correlates  with  KLK3  ( R 2  = 
0.6213  ,p  <  5  x  KT11)  and  KLK2  (R2  =  0.5893 ,p  <  5  x  1(T10), 
implying  that  KLK3e  may  coordinate  the  transcription  of  KLK3/2 
(Fig.  5B).  However,  we  do  not  completely  understand  how  KLK- 
associated  eRNAs  are  regulated  in  castration-resistant  prostate 
cancer  (CRPC),  although  it  is  well  accepted  that  the  AR-mediated 
transcriptional  program  remains  active.  This  observation  can  be 
explained  by  the  activity  of  other  nuclear  receptors,  because  the 
core  sequence  of  ARE  can  be  also  be  recognized  by  glucocorticoid 
receptor  (GR),  progesterone  receptor  (PR),  and  mineralocorticoid 
receptor  (MR)  (27).  More  importantly,  it  has  been  shown  that  AR 
inhibition  resulted  in  GR  up-regulation  in  a  subset  of  prostate 
cancer  cells  supporting  a  mechanism  of  bypassing  AR  blockade  via 
an  alternative  nuclear  receptor  in  CRPC  (28).  Thus,  it  is  possible 
that  GR,  PR,  or  MR  may  mediate  the  expression  of  KLK3e, 
KLK3,  and  KLK2,  attributing  to  adaptive  or  compensatory  effects 
generated  during  CRPC  development. 

In  view  of  the  regulatory  potential  of  IncRNAs  on  gene  regula¬ 
tion  in  trans  (29-31),  eRNAs  may  be  part  of  coactivator  complexes 
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and  regulatory  factors  that  participate  in  the  transcriptional  regu¬ 
latory  network.  Although  the  entire  circuitry  of  AR-dependent 
genes  is  unlikely  controlled  by  a  single  eRNA,  presumably,  there  are 
more  androgen-induced  eRNAs,  and  there  could  be  redundancy  in 
their  enhancer  activity.  It  is  conceivable  that  multiple  eRNAs  pro¬ 
duced  from  distinct  enhancer  elements  tether  with  their  own  tran¬ 
scriptional  modules  that  operate  simultaneously  to  initiate  an  ac¬ 
curate  and  a  rapid  transcriptional  program.  Therefore,  we  propose 
that  an  eRNA  may  function  as  a  scaffold  that  guides  an  AR- 
associated  protein  complex  to  target  chromatin  and  selectively 
enhances  DHT-stimulated  transcription  either  intrachromosomally 
(cis  activity)  or  interchromosomally  ( trans  activity).  This  positive 
regulatory  loop  may  provide  new  insight  into  RNA-dependent 
functional  effects  on  the  regulation  of  lineage-  or  tissue-specific 
gene  expression. 

Although  androgen  deprivation  therapy  (ADT)  is  designed  to 
block  the  androgen  signaling  pathway  and  is  the  standard 
treatment  for  advanced  prostate  cancer,  this  study  may  help  us  to 
understand  how  androgen  ADT  influences  the  activity  of  AR- 
bound  enhancers  through  altered  expression  of  eRNAs.  In  ad¬ 
dition,  the  differential  eRNA  expression  patterns  observed  may 
be  predictive  markers  of  ADT  efficacy  for  individual  prostate 
cancer  patients.  More  importantly,  eRNAs  may  represent  new 
targets  that  when  suppressed  may  improve  the  therapeutic  du¬ 
rability  of  ADT. 

Experimental  Procedures 

CHIP  Assay.  ChIP  experiments  were  performed  as  previously  described  (19) 
with  some  modifications.  Briefly,  10  million  cells  were  used  per  IP.  Cells  were 
fixed  with  1%  formaldehyde  solution.  DNA  was  sonicated  and  subjected  to 
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Experimental  Procedures. 
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